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Abstract 

To account for phenomenological theories and a set of invariants, stress and 
strain are usually decomposed into a pair of pressure and deviatoric stress 
' (-h ' 1 and a pair of volumetric strain and deviatoric strain. However, the conven- 

tional decomposition method only focuses on individual stress and strain, 
so that cannot be directly applied to either formulation in Finite Element 
Method (FEM) or Boundary Element Method (BEM). In this paper, a sim- 
pler, more general, and widely applicable decomposition is suggested. A new 
decomposition method adopts multiplying decomposition tensors or matrices 
to not only stress and strain but also constitutive and compliance relation. 
With this, we also show its practical usage on FEM and BEM in terms of 
tensors and matrices. 

Keywords: Pressure, hydrostatic pressure, Deviatoric stress, Volumetric 
strain, Mean strain, Deviatoric strain, Stress decomposition, Strain 
decomposition, Constitutive decomposition, Compliance decomposition, 
Multiplication decomposition, Decomposition multiplier 

1. Introduction 

In many references (see Fung (1965); Fung and Tong (2001); Gurtin (1981); 
Richards Jr (2000); Bower (2009)), stress and strain decomposition are usu- 
ally used for a set of invariants. From a hydrostatic stress tensor (or volu- 
metric stress tensor) p, the 1st stress invariant Ii is given by 
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h = 3p (1.1 

ress 

stress invariants are given by 



Also, from a deviatoric stress tensor s^, the second and third deviatoric 



■h — - Sij Sij (1.2) 

^3 = s y Sfcj (1-3) 

In 1.1-1.3, a hydrostatic stress tensor p and a deviatoric stress tensor sy 
result from the decomposition of stress tensor a if 



<y%i = Sij+pdij (1.4) 
P = -3- (1-5) 

Similarly, from a volumetric strain tensor 6m {^m = £kk/3) and a devia- 
toric strain tensor (e 4 - = ejj — e^/ <5jj ), we have strain invariants Jj, JJ, 
and J| as 



/i e = 3e M (1-6) 
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^3 = o e *j e jfe ( 1 - 8 ) 



The other general use of stress and strain decomposition can be found in 
phenomenological theories (e.g., Houlsby and Puzrin (2002); Perzyna (1966); 
Dunne and Petrinic (2005); Simo and Hughes (1998); Lubliner (1990); Leigh 
(1968)). For example, von Mises yield function / that describes elasto- 
plasticity is written as 



/ = -h-k 2 (1.9) 

where, k 2 = \ > a\. Also, the flow potential ip in Perzyna formulation that 
describes rate-dependent plasticity is written as 
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where, t] is a viscosity and (•) represents the Macaulay bracket (ramp 
function) . 

From 1.9-1.10, the flow rules for plasticity and viscoplasticity yield 



" i j 



A 



df_ 

doij 



A s,; 



(1.11J 



and 



6ij d Sij 

dip d.J 2 
dJ 2 dsij 



2 V \ y/J 2 

In 1.11-1.12, e?-, A, and e*J represent the plastic strain rate, plastic mul- 
tiplier, and viscoplastic strain rate, respectively. 

Despite such fundamental uses of decomposition in mechanics, the con- 
ventional decomposition method is not directly adopted to the formulation 
in FEM and BEM due to the lack of the generalized decomposition method 
including decomposition of constitutive or compliance relation. As we shall 
see, through the multiplication decomposition that developed here, one can 
directly decompose not only stress and strain, but also constitutive relation 
and compliance relation. The method is so simple and general that just a 
decomposition multiplier is applied to all cases mentioned here. And as an 
example, we will also show its applications to the formulation in FEM and 
BEM for elastostatics. 



2. Decomposition in tensor forms 

In this Section, we show the multiplying decomposition method for stress, 
strain, constitutive relation, and compliance relation in tensor forms. For 
canonical example, we also show its application to linear isotropic elasticity. 



2.1 Decomposition of stress and constitutive relation in tensor form 

2.1. Decomposition of stress and constitutive relation in tensor form 

Conventional stress decompositions such as 1.4-1.5 can be written in 
sor form as 



a 



ij 



fiki^lj — -&lk&ij ) °~kl 



and 



In 2.1-2.2, if we let 



and 



M fjki = SkiSij - -SikSij 



M ijkl = T^lkSij 



then, 2.1 and 2.2 can be written as: 

Sij = Mf jkl o k i 



and 

Pij = M ijki^ki 

In 2.6, represents 

Pij = pSij 

Such decompositions as 2.5 and 2.6 can be easily checked through 



2.1 Decomposition of stress and constitutive relation in tensor form 
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&ij &ij ~r Pij 

M? jkl a kl + M? jkl a kl 
= {5ki5ij)<Jki 

= (Tn (2.8) 
and, by them, we can also decompose the constitutive relation Cijki in 

&ij = Cijkl€kl (2.9) 
That is, by 2.9 and 2.5, the deviatoric stress can be written by 

Sij = M? jM a M 

= ^ijkl {C 'klmn^mn) 

= (Mfj k iCklmn) tmn (2.10) 

and, by 2.9 and 2.6, the pressure py can be written by 

Pij = M ijki^ki 

= ^ijkl {Cklmn^mn) 

= {^ijkl^klmn) tmn (2-11) 

Then, noting Mf jkl C klmn and M v ijk[ C klmn as 

Cjjmn = Mfj kl C k i mn (2.12) 

and 

Cij mn = Mi jkl Cklmn (2-13) 



yields stress decompositions 



Sij — C i j mn e rnn (2-14) 
Pij = C^ mn e mn (2-15) 



2.2 Decomposition of strain and compliance relation in tensor form 6 



in terms of decomposed constitutive tensors Cfj mn and C^ mn . 

Such constitutive decompositions as 2.12 and 2.13 are also checked through 



Cijmn + ^ijmn ~ (^ijkl + ^ijkl) Cklmn 
= 5ki5ljCklmn 

Cijmn (2.16) 

2.2. Decomposition of strain and compliance relation in tensor form 
Conventional strain decompositions are given by 

1 



4 = 1^ 



By following similar decomposition procedures to stress, the multiplying 
decomposition can define strain decompositions as 



4 = M: jkl e kl (2.19) 

e% = Mf jkl e kl (2.20) 

and, by them, we also can decompose the compliance relation D^i in 

e ij = Dijkio~ki (2-21) 

into 

D ijmn — M ijkl D klmn (2.22) 

and 

D ijmn = M ijkl D klmn (2.23) 



In 2.20, represents 



e# = e M 5 tJ (2.24) 



2.3 Example for linear isotropic elasticity 
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and equations 2.19-2.20 can also be written by 



e ij ~ ^ijmn a mn (2.25) 
e ij = ^ijmn^ran (2.26) 

in terms of the decomposed compliance tensors E)f- mn and D^ mn . 

Such compliance decompositions as 2.22 and 2.23 are also checked through 



Dijmn + ^ijmn ~ {^ijkl + ^ijkl) Dklmn 

Dijmn (2.27) 

2. 3. Example for linear isotropic elasticity 

For linear isotropic elasticity, the constitutive relation Cijki is given by 

Cijki = X5 k i5ij + 2/J,5 ki 5tj (2.28) 

By the suggested decomposition method, the constitutive decomposition 
tensors C^ mn and C^ mn are expressed as 



Cijmn — Mfj kl Cklmn 

8ki<5lj — ^Ik^ijj (^mn^kl + ^^mk^nl) 
2 

= (38 mi 6 n j - SijS mn ) (2.29) 

and 



Cijmn — M^ kl Cklmn 
1 



-(3A + 2n)6ij6 inn 

A + -ft) 5ijS mn 
KSij5 mn (2.30) 



2.3 Example for linear isotropic elasticity 
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where /x and K are shear modulus and bulk modulus, respectively. 

Equations 2.29-2.30 can also be checked through the conventional de- 
composition method. That is, from the elastic constitutive relation 2.28, we 
have 



= (A5 mn 5jj + 2fiS m iS n j) € mn (2.31) 
Letting ij — > kk in 2.31 yields 

o'kk = (3A + 2/x)e fcfc 

= 3Ke kk (2.32) 

Then, pjj that given by 2.7 is expressed as 

Pij — U 

= Ke kk 5ij (2.33) 
Also, given by 2.1 is expressed as 

&ij Pij 

&kk r 

= (^kiSij + 2fiS k iSij) e H - - ((3A + 2u)e kk ) 8^ 
2 

= (3eij - 8ije kk ) (2.34) 

Thus, from 2.33 and 2.34, we can verify the constitutive relations such as 

Pij = C i j mn e mn (2.35) 

and 

Sij = Cfj mn e mn (2.36) 

Similarly, the compliance tensor D^ k i for linear isotropic elasticity, which 
is given by 
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Dijki = — ((1 + v)5 k i5ij - v5 kl 5ij) 



(2.37) 



can be decomposed into 



D 



tjmn 



M ijkl D klmn 

2^ i^im^ jn T^ij&mn 



(2.38) 



and 



D v - 

i j ran 



3K 



[2.39) 



Such equations as 2.38 and 2.39 can also be checked through the conven- 
tional decomposition method. 

3. Decomposition in matrix form 

In this Section, we show the multiplying decomposition method for stress, 
strain, constitutive relation, and compliance relation in matrix forms. 

3.1. Decomposition of stress and constitutive relation in matrix form 

With Voigt notation, the decomposition multiplier 2.3 can be expressed 
in a matrix form as 
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(3.1) 



Also, the decomposition multiplier 2.4 can be expressed in a matrix form 

as 



3.1 Decomposition of stress and constitutive relation in matrix form 10 



1 

3 



1110 

1110 

1110 









(3.2) 



Then, we can express the deviatoric stress in 2.5 and the pressure in 2.6 
in a vector form s and p as 



and 



s=M d cr 



p=M"cr 



(3.3) 



(3.4) 



where, the Cauchy stress cr can also be expressed in terms of the consti- 
tutive relation matrix C as 



(T=Ce 

Substituting 3.5 into 3.3 and 3.4 yields 



(3.5) 



s = C d e (3.6) 
p = C v e (3.7) 

where, the constitutive relation C is decomposed into the matrices C d 
and C v : 



C d = M d C (3.8) 
C v = M v C (3.9) 

Matrix decompositions 3.3-3.5 can be checked through 



cr = s-\-p 

= M d (T + M v a 

= (M d +M v ) cr 

= Icr 



= cr 



(3.10) 



3.2 Decomposition of strain and compliance relation in matrix form 11 



where, I in 3.10 is the identity matrix of size 6. 

The matrix decompositions 3.6-3.9 can be also checked through 



cr = s+p 

= C d e + C v e 

= (M d C+M v C) e 

= (M d + M v ) Ce 

= Ce (3.11) 

3.2. Decomposition of strain and compliance relation in matrix form 

Similarly, we have the matrix decompositions for strain and compliance 
relation as 



M d e 
M v e 
D d cr 
D V (T 



(3.12) 
(3.13) 
(3.14) 
(3.15) 



and 



D d = M d D (3.16) 
D v = M v D (3.17) 

In 3.12-3.17, e d ,e v , D d , and D v represent the deviatoric strain, the vol- 
umetric strain, the deviatoric compliance relation, and the volumetric com- 
pliance relation in a vector and matrix form, respectively. 



4. Properties of multiplying decomposition method 

So far, we have shown the multiplying decomposition method for stress, 
strain, constitutive relation, and compliance relation in both tensor and ma- 
trix forms. The method just takes the multiplication decomposers Mfj kl (or 
M d ) and M^ kl (or M v ) and we will see their properties in this Section. As we 
shall see, the multiplication decomposition is also valid to decompose strain 
energy density and compatible with physical meaning. 



4.1 Decomposition of strain energy 
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4-1. Decomposition of strain energy 

With the multiplying decomposition method, strain energy density u can 
be written as 



u — - CijM €ij Cm 

= - {sM+Pki) (4 + 4f) 

= 2 Skie ' kl+ 2 Pkie * 1 ^ 
While deriving 4.1, we use the relations 



sw<# = (4.2) 
Pw4 = (4-3) 



since we have 



s ki e kl 



Ski tM&kl 





(4.4) 



and 



= V ( e kk — tM&kk) 

= (4.5) 

In 4.1, the multiplying decomposition method also decomposes the strain 
energy density u into deviatoric strain energy density s ki e ki) an d volu- 
metric strain energy density (| p^ ej^J. Such strain energy density decom- 
position can also be derived from u — | cr^ u^. 



4.2 Compatibility with physical meaning 
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4-2. Compatibility with physical meaning 

With the multiplying decomposition method, one may consider this can 
be also used for 

• the evaluation of total stress from either deviatoric stress (strain) or 
pressure (volumetric strain) 

• the evaluation of total strain from either deviatoric strain (stress) or 
volumetric strain (pressure) 

However, this is not possible because both the decomposition matrices M d 
and M v are singular. The eigenvalues of M d are Ai = and A2 = 1 with 
multiplicity 5 while the eigenvalues of M v aie Ai = with multiplicity 5 
and A2 = 1. In fact, the decomposition matrices M d and M v are projec- 
tion matrices. And having the eigenvalue of multiplicity 5 in M d and M v 
are natural since M d leaves only five deviatoric stress (or deviatoric strain) 
components and M v leaves only one pressure (or volumetric strain) compo- 
nent when multiplied to stress (or strain). Thus, we cannot have [iVf d ] 
and in 3.3-3.4 and 3.12-3.13. Consequently, we also cannot have 

[C d ]~\[CT\ l Dd ]~\ and [D v ]' X \n 3.8-3.9 and 3.16-3.17. 

In physical viewpoint, the properties of multiplying decomposition can 
be interpreted as: 

• Total stress results in both deviatoric and volumetric strain but not 
vice versa 

• Total strain results in both deviatoric stress and pressure but not vice 
versa 

5. Applications 

The multiplying decomposition method can be directly used to formulate 
both FEM and BEM. In this Section, we show its canonical application to 
elastostatics. 

5.1. Finite element formulation 

In FEM for elastostatics (see Strang and Fix (1973); Cook et al. (2002); 
Bathe (1996); Hughes (2000); Braess (2001)), the element stiffness matrix 
K, is determined by 



5.1 Finite element formulation 
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K e = [ B T C B dVt (5.1) 

where, B,C, and Q e are the strain-displacement matrix, the constitutive 
relation matrix, and the domain of the element, respectively. 

After assembling all the stiffness matrices of each element and accounting 
for boundary conditions to describe the given structure, the global system of 
equations is given by 

K u=f (5.2) 

In 5.2, K,u and / represent the global stiffness matrix, nodal displace- 
ments and nodal forces in a vector form, respectively. 

With the adoption of the multiplying decomposition, we can reformulate 
the element stiffness matrix K e in 5.1 as 

K e = [ B T (C d +C v )Bdn 

= j B T C d BdQ+ I B T C V B dVt (5.3) 
= K d + K v 

e e 

where K d and K v e are decomposed element deviatoric stiffness matrix 
and element volumetric stiffness matrix. 

Also, the global system of equations can be written as 

(K d + K v ) u=f (5.4) 



K K' 1 K' (5.5) 

As in 5.5, with the multiplying decomposition, the global stiffness matrix 
K is decomposed into the global deviatoric stiffness matrix K d and the 
global volumetric stiffness matrix K v , where K d results from assembling all 
the K d and K v results from assembling all the K v . 

Separate evaluation of K d and K v by the multiplying decomposition 
allows "selective integration technique" in a direct way and some useful ap- 
plications are listed below. 



5.2 Boundary element formulation 
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• Using the reduced integration for volumetric part and full integration 
for deviatoric part: in (nearly) incompressible cases, selectively reduced- 
integration technique (SRI) is used for 1st order linear isoparametric 
elements (4-node elements in two dimensions and 8-node elements in 
three dimensions), to prevent mesh locking and get accurate solution for 
(nearly) incompressible cases (see Malkus and Hughes (1978); Hughes 
(2005); Doll et al. (2000); Liu et al. (1994, 1998)). 

• Using the reduced integration for deviatoric part and full integration 
for volumetric part: shear locking can be also prevented by SRI. (See 
Bathe (1996); Cook et al. (2002); Braess (2001); Hughes (2000)) 

• Formulation of plasticity: when von Mises plasticity occurs (see 1.9, 
1.11), K d is evaluated iteratively while K v remains constant since 
e kk = (See Dunne and Petrinic (2005); Lubliner (1990)). With the 
adoption of the multiplying decomposition, we can separately evaluate 
K d without evaluating K v during plastic evolution (computation can 
be reduced). 

5.2. Boundary element formulation 

In BEM, an integral equation for interior stress in elasticity can be written 
as (see Banerjee (1994); Beer et al. (2010); Ang (2007); Wrobel and Aliabadi 
(2002)) 

/ {C^(x,t)t k (x)-Fk(x,t)u k {x))dS(x) (5.6) 
Js 

where G^- and F ki - are traction and displacement kernel function for 
stress respectively, and S means the boundary of a given problem domain. 
With the multiplication decomposition, deviatoric stress can be derived 

as 

= Mf jmn f (Cfy(x,t)t k (x)-Fh{x,t)u k {x))dS{x) (5.7) 
Js 

Similarly, pressure can be derived as 
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Pmn {i) = M v ijmn f {(^(x,$)t k {x)-F^{x^)u k (x))dS(x) (5.8) 
Js 

Similarly, decomposed strains are as follows 

e' mn (Z) = Mf 3mn I (G e Hj (x, £)t k {x) - F^(x, £)u k (x)) dS(x) (5.9) 
Js 

e%n(t) = M? jmn [ (Gl i:j (x^)t k (x)-F^(x^)u k (x))dS(x) (5.10) 
Js 

where and F^- are traction and displacement kernel function for 
strain respectively. 

6. Conclusions 

A simple, clear, and widely applicable way to decompose stress/strain 
and constitutive/compliance relations is suggested in both tensor and ma- 
trix forms: multiplication decomposition. The method is also applicable to 
decompose strain energy density along with proper physical meaning. 

We consider here the application of multiplying decomposition to elasto- 
statics in FEM and BEM formulation, which illustrates the elegance of this 
approach. Clearly, however, the multiplying decomposition is quite general 
and can be applied readily to elastoplasticity, viscoplasticity, fluid mechanics 
and more broadly throughout mechanics. In addition, we anticipate that 
the multiplying decomposition method developed here will provide an inter- 
esting foundation for the development of novel analytic and computational 
methods. 
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